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Schramm-Loewner Evolution (SLE) is a stochastic process that helps classify critical statistical 
models using one real parameter k. Numerical study of SLE often involves curves that start and 
end on the real axis. To reduce numerical errors in studying the critical curves which start from the 
real axis and end on it, we have used hydrodynamically normalized SLE(fi;, p) which is a stochastic 
differential equation that is hypothesized to govern such curves. In this paper we directly verify this 
hypothesis and numerically apply this formalism to the domain wall curves of the Abelian Sandpile 
Model (ASM) (k = 2) and critical percolation (k = 6). We observe that this method is more reliable 
for analyzing interface loops. 
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I. INTRODUCTION 

Classification of probability measures on random 
curves in simply connected domains of the complex plane 
was first proposed by O. Schramnii. According to this 
idea one can describe the interfaces of two dimensional 
statistical models via growth processes named as SLE. 
This technique has become a powerful tool to study the 
macroscopic interfaces of two dimensional systems. An- 
other powerful tool for classification of 2d critical phe- 
nomena is conformal field theory (CFT). Hence a con- 
nection between the two methods is natural. The rela- 
tion between CFT and SLE was found in the pioneering 
works of M. Bauer and D. Bernard^,. Parallel to the ana- 
lytical achievements, a large amount of numerical works 
have been done using SLE to study the statistics of the 
critical interfaces of various statistical models. Since the 
interfaces of the critical models resulting from simula- 
tion or experiment are commonly loops, the customary 
method that is used to make curves going from origin to 
the infinity is to cut the loops from its middle and send 
the end point of the cut curves to infinity (infinity map- 
ping). This method is widely used e.g. in Turbulence^, 
ASM avalanche frontier— iso height lines of KPZ^ and 
W03^ , Ising modeli etc. But there is an important 
source of numerical error due to the infinity map which 
sends the end point to the infinity. This map enlarges 
the lattice constant in the regions close to the end point 
of the curve and this causes numerical large errors. To 
avoid these errors, we should develop a method that is 
free of infinity mapping. 

The new method we introduce in this paper, is based 
on a variant of SLE^, known as SLE(«;,/9). SLE(K,yo) is 
a generalization of the SLE^ that describes the curves 
that are self similar but are not conformally invariant 
i.e. there are some preferred points on the domain that 
the SLE curve is growing. It has been claimed^ that 
the hydrodynamically normalized stochastic differential 
equation governing the mappings of the critical curves 
that start from the real axis and also end on it, is the 
SLE(ft, p) where p = k— 6. Such maps are the ones we use 
to develop our numerical method. Therefore we not only 
introduce a new more reliable numerical method to ana- 
lyze SLE curves, but also we will have a numerical check 



for the application of SLE(k, p) formalism to the curves 
going from the real axis to the real axis as stated above. 
We show that this formalism is properly applicable and 
obtain more precise results for the ASM and percolation. 
For this end, we use the slit uniformizing map to obtain 
the driving function and use the Maximum Likelihood 
Estimation (MLE) to fit properly the resulting process to 
the Brownian Motion to find the parameters p and n and 
their stochastic errors. In section |TT] we briefly introduce 
chordal SLE. Sections IIIII is a brief introduction to the 
critical percolation and ASM respectively and SLE(k, p). 
In section ITVl we present our numerical results. 



II. SLE 

Critical behaviour of the two dimensional statistical 
models can be described by their geometrical features. 
In fact instead of studying the local observables, we can 
focus on the interfaces of two dimensional models. These 
domainwalls are some non-intersecting curves which di- 
rectly reflect the status of the system in question and sup- 
posed to have two properties: conformal invariance and 
the domain Markov propertyi^. Schramm- Loewner Evo- 
lution is the candidate to analyze these random curves 
by classifying them to the one-parameter classes (SLE^). 



A. Chordal SLE 

Let us denote the upper half plane by H and 7f as 
the SLE trace i.e. = {z <E H : Tz < t} and the 
hull Kt = {z G H : Tz < t} . SLE^ is a growth pro- 
cesses defined via conformal maps which are solutions 
of Loewner's equation: 



dt9t{z) 



9tiz)~^t 



(1) 



Where the initial condition is gt{z) = z and = \ffiBt 
is a real valued smooth function. For fixed z, gt{z) 
is well-defined up to time for which gt{z) = S,t- 
The complement Ht := H\Kt is simply connected and 
gt{z) is the unique conformal mapping Ht ^ H with 
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gt{z) = z + ^ + O(p-) as z — 7> CX3 that is known as hy- 
drodynamical normalization. One can retrieve the SLE 
trace by 7t = Unijj^o <?t~^(Ct + There are phases for 
these curves, 2 < k < 4 the trace is non-self- intersecting 
and it does not hit the real axis; in this case the hull 
and the trace are identical: Kt = 7t. This is called "di- 
lute phase". For 4 < k < 8, the trace touches itself and 
the real axis so that a typical point is surely swallowed 
as i — > oo and Kt ^ 74. This phase is called "dense 
phase" . There is a connection between the two phases: 
for 4 < K < 8 the frontier of Kt^ i.e. the boundary of 
Ht minus any portions of the real axis, is a simple curve 
which is locally a SLE;^ curve with k = — , that is it is 
in dilute phase. The main question "what is the relation 
between SLE and CFT" is answered by M. Bauer and 
D. Bernard^. They showed that the boundary condition 
changing (bcc) operator in SLE correspond to a degen- 
erate field with a vanishing descendant at level two and 
conformal weight h\.^2 ~ in CFT with central charge 
c„ = (&~K:)j^j^-'&) ^ tjjaj; |;j^jg relation is consistent 

with the duality expressed above: c(k) = c(k ~ 16/k). 

B. SLE(k,p) 

Now we will give a brief introduction to SLE(k, p). We 
define SLE(k, p) in the upper half plane. The parameter 
K identifies the local properties of the model in hand and 
corresponds directly to its central charge, the parameter 
p has to do with while the boundary conditions (be) im- 
posed. One be occurs at the origin (the starting point of 
the curve) and the other at some point x^o on the real 
axis. The stochastic equation governing on such curves 
is the same as formula ([!]) but the driving function has 
different form, namely; 

dit = V^dBt + -dt (2) 

Now consider a curve that starts from origin and end on 
a point on real axis (xoc)- Then by using the map (/) = 
Xaczjixoa — z\ One cau send the end point of the curve to 
the infinity. In this respect, the function ht = 4>ogtO(j)~^ 
describes chordal SLE. It is easy to show that the equa- 
tion governing on ht is dtgt = 2/{{(f)'{gt){(j){gt) -6)})- 
But it is explicit that this function is not hydrodynam- 
ically normalized. It has been shownH that if one uses 
another mapping gt = vtohtou~^ where u = 4>~^ and vt 
is any linear fractional transformation that make the cor- 
responding map hydrodynamically normalized, then the 
stochastic equation of gt is the same as Eq. ([T|) . In fact 
this procedure leaves the Eq. ([T|) unchanged but leads 
the driving function to have a drift terroi^: 

d^t = y/^dBt + —^^—^—dt (3) 

In the other words, this stochastic function is the driv- 
ing function of the SLE(k,p) with p = k — 6. Thus for 



the critical curves from boundary to boundary, the cor- 
responding driving function acquires a drift term. This 
generalization of SLE can be generalized further to have 
multiple preferred real axis points. For review see [13]. 



III. STATISTICAL MODELS 

In this section we briefly introduce some statistical 
models to be simulated in the next section. To cover the 
argument properly we choose models from two regimes 
of K that are ASM with k = 2 < 4 (dilute phase) and 
percolation with 4 < k = 6 < 8 (dense phase). In below 
we give a short review of the two models. 



A. Critical Percolation 

Let (H, a, b) be an hexagonal lattice domain on the 
upper half plane with the boundary condition changes 
on two points on the real axis namely a and b i.e. color 
the hexagons on the boundary in black and white so that 
the resulting boundaries have two boundary condition 
changes on points a and b as indicated in the Fig. [T] . A 
configuration is a choice of color for inner hexagons which 
becomes black or white with the probability p and 1 — p 
respectively. One can easily identify unique domainwalls 
separating white and black sites from each other. Each 
configuration defines an interface, i.e. the unique path 
from a to b in H such that the hexagon on the left (resp. 
right) of any of its edges is black (resp. white), see Fig. 
[1] Hence the probability distribution on configurations 
induces a probability distribution on paths from a to b 
in H. There is a critical probability p = pc = 0.5 in which 
the size of the resulting clusters is of order of the size 
of lattice and the clusters are self similar. We consider 
here the critical interfaces of percolation which is known 
that are SLE(k ~ 6). Percolation enjoys also an impor- 
tant property namely locality, that is, the evolution of 
the SLE curve is insensitive to the boundary conditions 
and the events that take place in the H \ Kt up to the 
time t. Due to the locality, the interfaces of the perco- 
lation can be generated by a simple growth process: in 
the n'th step, when the tip of the curve reaches to an 
inner hexagon, then using a fair coin, color this site in 
black or white. The upshot is that the n + I'th step of 
the resulting domainwall will be the edge of the hexagon 
whose adjacent faces have different colors. By this simple 
rule, we can generate domainwall samples that start from 
the origin and end on the arbitrary point on the lattice 
boundary. Locality implies that this growth process is 
insensitive to the fact that "where the curve will end". 
It is guaranteed in Eq. ^ which implies p = k — 6 = 0, 
saying that the drift term in this equation vanishes. So 
these stochastic curves evolve like the chordal one and 
are insensitive to the boundary conditions as expected. 
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FIG. 1: A typical configuration of tlie percolation starting 
and ending on the real axis. 



B. ASM 

Sandpile models have been introduced by Bak, Tang 
and Wiesenfeld^ as an example for a class of models 
that show Self-Organized Criticality(SOC). These mod- 
els show critical behaviours, without tuning external pa- 
rameters such as temperature. The abelian structure of 
this model was first discovered by D. Dhar and named as 
Abehan Sandpile Model (ASM)iS. Despite its simplic- 
ity, ASM has various and interesting features and many 
different analytical and numerical works have been done 
on this model. For example different height and clus- 
ter probabilities^, its connection with spanning trees^i, 
ghost modelsi^; q-state potts modeU^, etc. For a good 
review see [16]. 

Consider the ASM on a two dimensional square lattice 
L X L. To each site i, a height variable hi is assigned 
taking its values from the set 1,2,3,4 the number of sands 
in this site. The dynamics of this model is as follows; 
in each step, a grain is added to a random site i i.e. 
hi hi + 1: if the resulting height becomes more than 
4, the site toppcls and loses 4 sands, each of which is 
transfered to one of four neighbors of the original site. 
As a result, the neighboring sites may become unstable 
and topple and a chain of toppling may happen in the 
system. In the boundary sites, the toppling causes to one 
or two sands to leave the system. This process continues 
until the system reaches to a stable configuration. Now 
another random site is selected and the sand is released 
on this site and the process continues. The movement on 
the space of stable configuration lead the system to fall 
in a subset of sets of configurations after a finite steps, 
named as the "recurrent states" . For details see [16]. This 
model is related to Potts model with q = and CFT 
with c = 0. For a lattice with d neighboring sites, the 
toppeling occurs when hi > d, then the original site will 
lose d grains and the height of each of its neighbors will 
increase by 1. It is obvious that for a triangular lattice 
d = 6. 



FIG. 2: The schematic view of triangular lattice with two 
sublattices A and B. 



IV. NUMERICAL METHODS AND RESULTS 

To use Eq. ([ij for the statistical models on the lattice, 
one has to discretize this equation. For this purpose it 
is customary to assume an especial function form for the 
driving function in each discrete interval and find the cor- 
responding uniformizing map in that interval and using 
Gt„ ~ Gt^_ioGt^_20...oGti send every point on the -ft to 
the real axis step by step. If we assume that ^ is constant 
in an interval (^n), the corresponding uniformizing map 
is 



Gt„ (z) = + Viz - e„)2 + 4t„ (4) 

It is notable that gt{z) = z + ^ + 0{^) as z — ^ oo. 
The method that is widely used to apply this formalism 
to the interface loops for the statistical models on lattice 
is as follows: one cuts the loops horizontally from its 
middle to generate curves starting from the origin and 
ending at some point on the real axis. Then by mapping 
the end point of the curve to infinity and applying the 
chordal SLE i.e. Eq. ([1]) one can obtain the underlying 
driving function. For the simulation, we consider the 
triangular lattice and its dual lattice (honeycomb lattice) 
and also consider two sublattices A and B as indicated 
in Fig. [2] Thus the position of each point can be coded 
in (n,m,A/B) where n and m shows the position of the 
lattice point and A or B shows its basis. It is obvious that 
on this lattice, each point has 6 neighbours. Consider 
that we have ensembles of interface loops of some model, 
then by the process described above we have ensembles 
of the curves that start and end on the real axis. Thus 
when we have the set of discrete points {zi, Z2, z^, Zn\, 
the first unformizing map is of the form of Eq. Q with 
the parameters 5ti = |;(Im[zi])^ and = Re[zi] and 
after applying this mapping we continue this process for 
Z2. The process continues until every point of the curve 
is sent to the real axis. To analyze the resulting driving 
function and use Eq. (j?]), we discretize this equation. 

<5Cn = \f^5Bn + 7 z5tn (5) 



FIG. 3: A sample graph of the percolation going from real 
a:xis to itself. 



FIG. 4: ^ graph vs k, it has a minimum at «; = 6.08 and 
po = -0.08. 



For having a pure Brownian motion, wc should write the 
Eq. (6) of the form; 



(6) 



Where p = k—6. After doing this process, we can analyze 
the resulting quantity on the right hand side of Eq. (|6]) 
to see if this stochastic process is Brownian Motion. The 
best fit to the Brownian Motion, gives us the amount of 
K. For this, we used the Maximum Likelihood Estima- 
tion (MLE) method. In this method, the best value for a 
model parameters (here k and po = —p) are which min- 
imize the function and their precisions can be obtained 
from the distribution function of the function by cal- 
culating the ratio of the area under the resulting diagram 
to the whole area under it. 




A. Critical Percolation 

Consider the interfaces of critical percolation. A typ- 
ical critical percolation curve is shown in Fig. [3] that 
starts from origin and ends at the point 5200 on the real 
axis. Due to its locality property, percolation is so hard 
to be controlled to end on a certain point, in this figure 
it reaches the real axis after 850000 steps. We simulated 
over 10000 curves with the typical length 50000. After 
simulating such curves we obtained as described above. 
The Fig. shows the behaviour and the minimum of the 
function versus k which is obtained by fitting {Bf) 
versus t. The global minimum of this curve is in the 
point K = 6.08 and po = 6 — k = —0.08 in agreement 
with the expected values (k = 6 and po = 0). This result 
justifies the locality property of the percolation as stated 
above. 

Using the distribution of we obtain the precision 
of the resulting parameters that is, with the probability 
60% K 6.08 ± 0.07 and po = -0.08 ± 0.07. The plot of 
{B'I) versus t graph for the obtained k, and po is indicated 
in Fig. [5] We see that this graph is linear with the slope 



FIG. 5: The graph of the obtained process {B'i) versus t for 
percolation. 



0.98 ± 0.02. The inner graph shows the distribution of 
BtS for t ~ 0.625 that is properly fitted to the Gaussian 
function with a ~ 0.625 this ensures us that this process 
is also a Brownian Motion. 



B. ASM 

Consider now avalanche frontiers of ASM. By adding a 
grain to the saturated sandpile and making it unstable, 
avalanches occur and their frontier form a loops with dis- 
crete points. Fig. |6]shows a typical interface of the ASM 
that goes from real axis to the real axis. As stated above 
from the theory one expects that po = 6 — k = 4. To test 
this, we let po to be a free parameter in this case. The re- 
sulting amounts for the ASM on the 2048 x 2048 triangu- 
lar lattice over 10000 interface curves are k = 1.95 ±0.05 
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FIG. 6: A sample interface of ASM starting from origin to 
the real axis. 
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that is the famous property of the Brownian motion. In 
the inner graph of Fig. [7] the distribution of Bt at time 
t = 0.125 and also the Gaussian distribution exp[— -^J 
with (7 = 0.125 are indicated in the same time. 

As is seen in the graph they do suitably fit to each 
other showing that the distribution of i?t's arc Gaussian. 
These evidences and the fact that {Bt) ~ show that 
this stochastic process is Brownian Motion and the driv- 
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FIG. 8: The contour graph for the (for fitting (Bt) with 
t) as a function of k and po. 



ing function obtained from hydrodynamically normalized 
uniformizing map process for this process has the form 
ofEq. (ED. 



FIG. 7: The {B^) graph vs t that is shown is linear with a good 
slope 1.00 =F 0.05 for ASM. The inner graph is the distribution 
of Bt for t = 0.125. The red graph is the numerical result 
and the blue one is the Gaussian distribution exp[— -^J with 
(7 = 0.125. 

and Po = 3.5 ± 0.5 with the probability 60%. Fig. E 
indicates the the contour plot of in terms of k and pa 
in which the probability of finding parameter values in 
various regions are indicated. In Fig. [7] we have shown 
the graph (B^) versus t for the parameters considered 
above {k = 1.95 and po = 3.80). It is seen that the graph 
is properly linear with the slope 0.98 ± 0.05 {{Bf) ~ t) 



V. CONCLUSION 

In this paper, we analysed the statistics of the curves 
that go from real axis to itself using the formalism of hy- 
drodynamically normalized SLE(k,/9) which is expected 
to describe them on the upper half plane, to test its valid- 
ity and make the numerical calculations more exact. For 
this end we considered avalanche frontiers of the criti- 
cal ASM and the critical percolation. Then fitting to the 
Brownian Motion and using the Maximum Likelihood Es- 
timation, we calculated the best values of parameters i.e. 
corresponding k and po = —p and their precisions. In 
this way wc introduced a new method to analyze such a 
curves that is more precise and more reliable. 
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